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Testing Theoretical Models for the Higher-Order Moments 
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ABSTRACT 

Using high-resolution N-body simulations, we test two theoretical models, based ei- 
ther on spherical or on ellipsoidal collapse model, for the higher-order moments of the 
dark matter halo distribution in CDM models. We find that a theoretical model based 
on spherical collapse describes accurately the simulated counts-in-cells moments for 
haloes of several mass ranges. It appears that the model using ellipsoidal collapse 
instead of spherical collapse in defining dark haloes is unable to improve the models 
for the higher-order moments of halo distribution, for haloes much smaller than M* 
(the mass scale on which the fluctuation of the density field has a rms about I). Both 
models are particularly accurate for the descendants of haloes selected at high red- 
shift, and so are quite useful in interpreting the high-order moments of galaxies. As 
an application we use the theoretical model to predict the higher-order moments of 
the Lyman break galaxies observed at z rs 3 and their descendants at lower redshifts. 

Key words: Galaxies: formation - galaxies: clustering - galaxies: haloes - cosmology: 
theory - dark matter 



1 INTRODUCTION 

In the standard scenario of galaxy formation, it is assumed 
that galaxies form by t he cooling and condensation of gas 



199e ; LJing 1998 



poib 



199£ ; polberg et al 



S heth 

pooop . 



Tormen 199E; 



Governato et al 



withi n dark haloes (e.g. White fc Rccs 



1978 



White & Frenk 



|l99l| ). The problem of galaxy clustering in space can then 
be approached by understanding the spatial distribution of 
dark haloes and galaxy formation in individual dark haloes. 
This approach is very useful for the following two reasons: (i) 
the formation and clustering properties of dark haloes can 
be modelled relatively reliably because of the simple physics 
involved (gravity only), (ii) realistic models of galaxy for- 
mation in dark haloes ca n now be constructed using either 



The results show that the model 
proposed by Mo & White works reasonably well over a large 
range of halo masses. However, significant discrepancy be- 
tween model and simulations results was found for low-mass 



haloes ( |Jing||l998| ; |3hcth fc Tormerj|l999| ). gheth et alj (gOOjj ) 



semi-analytic models (e.g. Kauffmann et al. 1999; Cole et al 



200C ; Somerville & Primack 



lations (e.g. Benson et al. 2001). Indeed, attempts have been 
made to use theoretical models of halo clus tering to under 



1999) or hydrodynamical simu- 



Jing et al. 1998 



.; M 

Peacock fc Smith||2000| |Scljak| pOOOp . Most of these invest! 



stand clustering properties of galaxies (e.g. Mo et al 
~ ~Fryl |2000|; fj coccimarro et al 



1997 



2001 



gat ion s use the theoretic al mod els presented in 



(1996) and in Mo et al. (1997) (hereafter MJW) to calcu 



Mo & White 



late the second-order and higher-order correlations of dark 
haloes. These models are based on the Press- Schechter for- 



(hereafter SMT) suggested that the discrepancy at the low- 
mass end may be due to the fact that the model consid- 
ered by Mo and White assumes spherical collapse for the 
halo formation while the collapse in a realistic cosmological 
density field may be better approximated by an ellipsoidal 
model. Indeed, SMT found that, if an ellipsoidal model is 
used, better agreement between the model and simulations 
results can be achieved in both the halo mass function and 
the two-point correlation function for low-mass haloes. 

The performance of the MJW model for the higher- 
order moments of the halo distribution has been tested in 
their original paper using scale-free N-body simulations with 
relatively low resolution. Although their results show that 
the theoretical model matches the simulations results, the 
limited dynamical range in the simulations used by them 
does not allow one to test the model for a large range of 
halo masses. Furthermore, although the MJW model has 



& Cole 



1994) 



malism ( Pres s & Schechter 1974) and its extensions (Lacey 



The model prediction for the second moment, or the 
two-point correlation function, has been tested qu i te exten- 
sively by numerical simulations ( Mo fc White 1996 : Mo et al 



been extended to include ellipsoidal dynamics (Sheth et al 
2001), this extension has not yet been tested by simulations 
results. 

In this paper we use two sets of high-resolution simu- 
lations to test the MJW model and its extension. One set 
has a very large simulation box (and thus low mass resolu- 
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tion) which is used to control the finite- volume effect usually 
found in the analysis of higher-or der moments of the galaxy 
distribution (Colombi et al. 1994). The other set has smaller 
simulation boxes but much higher mass resolutions which al- 
lows us to test the models for low-mass haloes. 

The paper is organized as follows: The procedure to 
obtain the high order moments from counts-in-Cells is pre- 
sented in section 2 along with theoretical models of these 
moments for dark haloes. Analysis of the simulations data 
and the comparison of the theory with the simulations re- 
sults are presented in section 3. An application of the the- 
oretical models to the higher-order moments of the Lyman- 
break galaxies (L BGs) is done in section 4. Because LBGs 
are highly biased ( Adelberger et al. 1998 ), the skewness and 
kurtosis coefficients (S3 and S4) of LBGs and their descen- 
dants are significantly lower than those of the mass. 

Finally, section 5 contains a summary of our results. 



2 HIGHER-ORDER MOMENTS OF 
COUNTS-IN-CELLS 

2.1 Definitions 

The calculation of the counts-in-cells moments of a particle 
distribution and the relation of such moments to the corre- 
sponding moments of the underly i ng con tinuous density field 
are described in detail in 



Peebles (1980). We summarize the 



relevant formulae in the following. 

The j th central moment of counts in cells (we will use 
spherical cells whose radius will be denoted by R) of a point 
distribution is defined as 



1 M 



(i) 



where Ni is the number of particles counted in the i* 
sphere (cell), iV is the mean number of counts: N(R) = 
jj 12tLi Ni(R), and the summation is over the M sampling 
spheres. Notice that N is obtained directly from the counts. 

The connected moments, /a, are defined through the 
central moments as 



fJ-2 = m 2 , 
A*3 = m 3 , 
[i4 — m,4, - 



3ml . 



(2) 
(3) 
(4) 



These relations are written up to the 4th order because to 
this order they are relevant to our subsequent discussion. For 
a point process, the shot noise also contributes to the quanti- 
ties fij . These contributions become significant for small ra- 
dius where the mean count N is small and should be properly 
subtracted. If the particle distribution is a Poisson sampling 
of the underlying density distribution, we can make the fol- 
lowing subtractions to get the corrected connected moments: 



fc 2 — fi2 — N, 



fJ-s - Sn2 + 2N, 



It4 — fl4 — 6/I3 + 11/12 — GN. 



(5) 



(6) 



(7) 



These quantities are related to the volume- averaged corre- 
lation functions by 



(8) 



where 



dn...dr 3 W(n)...W(r 3 ) fc(ri, ...,!•,■), (9) 



and W(r) is a top-hat spherical window with volume Vw- 



2.2 Theoretical Model for the Higher-Order 
Moments of Dark Halos 



Mo et al. (1997) (MJW) have developed an analytical 



model for the hierarchical correlation amplitudes [Sj t h(R) = 
^j,h/^2~h ] f° r 3 = 3,4,5 in the quasi- linear regime, where 
the subscript h stands for quantities of dark haloes. In this 
model the statistical distribution of dark haloes within the 
initial density field, which is assumed to be Gaussian, is de- 
termined by an extension of the Press-Schechter formalism 
and the modifications of the distribution due to gravitation- 
ally induced m otions are treated by means of a spherical 
collapse model (Mo & White 1996). The main results from 



this model, which are relevant for our analysis, are summa- 
rized as follows. 

For the skewness and kurtosis (j = 3 ,4), following the 
same notation as in MJW, one has: 



b-^Sa + Sca), 



6- 2 (S4 + 12c 2 5 , 3+4c3 + 12c2 



(10) 



(11) 



where 53 and 5*4 are the skewness and kurtosis of the under- 
lying mass density field, Ck = bh/b, b = b\, and the constants 
bk are the coefficients in the expansion of the bias relation: 



00 . 



(12) 



where 8^ is the overdensity of haloes smoothed in a given 
window and 8 is the corresponding overdensity of mass. The 
coefficients bk for a halo with mass Mi corresponding to a 
linear overdensity 81, which collapses at redshift Z\ = S 1 /8 C — 
1 (with the critical overdensity for spherical collapse being 
8c — 1.686), are given by: 



6i = l + 



8x 



2(1 



6(02 + a 3 



■ CL2) 



81 



K-3), 



(13) 



(14) 



+ 



2 4 



61 

-6^i+3 
Si 



+ 3(1 +2a 2 ) 



K 2 - 3) 



(15) 



where vi = 8i/a(Mi) [with er(Mi) being the rms of the den- 
sity fluctuation given by the density spectrum linearly ex- 
trapolated to the present time], 02 = — 1| and as = ||i are 
coefficients in the expansion of So(8), the relation between 
the real mass overdensity 8 and the corresponding quantity 
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obtained using linear theory (see equation A4 in MJW). The 
above bias coefficients are for the present-day descendants 
(at redshift zo = 0) of haloes identified at redshift z\. It 
can be easily extended to the case where 21 > zo > 0. In 
this case, we replace <5i by 8\D(zo) / D(0) (where D(z) is the 
linear growth rate evaluated at redshift z) while keeping v\ 
unchanged. 



If £ n ellipsoidal model is used to define collapsed haloes, 



the coeffic ients bk take the following forms (Scoccimarro 



et al. |pq)0lD 



61 = l + e 1 +E u 



62 = 2(1 + oa)(ei + Ei) + e 2 + E 2 , 



(16) 



(17) 



63 = 6(a 2 +a 3 ) (ei+£;i) + 3(l + 2a2) (e 2 + E 2 ) + e 3 + E 3 , (18) 
where 



ei = 



o 2 



£3 



Ei 

Ei 

E x 
E 3 
Ei 



61 ' 
-gr(av -3), 

ap2 I 2 4 a 2 . Q \ 

— r=-(a f — 6<w +3), 
5 i 



2g/fo 



5i 



(a^ 2 )!>' 
^ +2ei 



^ = 4(p 2 



6pa^ 2 

1) H j2 H 3e a , 

"1 



(19) 



(20) 



(21) 

(22) 
(23) 



and a = 0.707, p = 0.3. These formulae reduce to the origi- 
nal MJW model for a = 1 and p = 0. 

Inserting the expressions for bk in equations (^) and 
( pr[ ) and taking S3 and S4 in these equations as the skewness 
and kurtosis of the mass distribution measured directly from 
the N-Body simulations, we can calculate the skewness and 
kurtosis for the distribution of dark haloes as predicted by 
the MJW model and its SMT extension (i.e. the model with 
ellipsoidal collapse). 



(see below). We compare the results with those obtained 
from the GIF simulations to see how comparisons between 
models and simulations can be made for simulations with 
relatively small volume. For the GIF simulations, we focus 
on the rCDM and ACDM models. The parameters charac- 
terizing the simulation s are summarized in tab le 1. Further 
details can be found in 



et al 



(1998) 



Kauffmann et al 



(1999) and Jenkins 



For each simulation there are several output files corre- 
sponding to different evolutionary times (redshifts) and for 
each of these output times there is a halo catalog contain- 
ing information about haloes identified using the friends- 
of- friends group-finder algorithm with a linking length 0.2 
times the mean interparticle separation. Only haloes con- 
taining 10 or more particles are included in the halo cata- 
logues. In the case of the GIF simulations at redshift zero, 
which are used to compare the predictions of the models 
for haloes with low masses, unbound objects are excluded 
from the catalogues (see section 3.2 for details). The phys- 
ical quantities available from each of these halo catalogues 
are: the index of the most-bound particle in the halo, which 
corresponds to the position of the halo as well as the central 
'galaxy' within it; the virial radius (R V i r ), defined as the 
radius (from the central particle) within which the overden- 
sity of dark matter is 200 times the critical density; the virial 
mass (M V i r ), which is the mass (or, equivalently, the total 
number) of dark matter particles within the virial radius; 
the circular velocity [V c = (GMrir/Eri,) 1 ' 2 ]. 

We have also generated several catalogues of the 
present-day positions of the central objects corresponding 
to the most-bound particles in haloes identified at an earlier 
epoch. This catalogues might be interpreted as 'galaxy cat- 
alogues' if we assume that the positions of galaxies at the 
present epoch correspond to those of the central particles 
within virialized objects identified at high redshifts. This 
concept is related to the assumption in models of galaxy 
formation that galaxies form by the cooling and condensa - 
tion of gas within dark matter haloes (White & Frenk 1991). 
However, this interpretation does not take into account sub- 
sequent galaxy mergers. 

We apply the counts-in-cells analysis described in the 
last section to the mass distributions and halo catalogues. 
To do this, we place spheres in a regular mesh of 30 3 centers 
and count the number of objects at each center over a set of 
concentric spheres which allows us to compute the desired 
statistical quantities at different radii. 



3 SKEWNESS AND KURTOSIS OF DARK 
HALOES IN N-BODY SIMULATIONS 

3.1 Simulations 

In the present analysis we use two sets of cosmological N- 
body simulations, which have been obtained as part of the 



3.2 Results 



VIRGO (Jenkins et al. 1998) and the GIF (Kauffmann et al 



199S) projects. These two sets of simulations differ in the 



size of the simulation boxes and in the mass resolution, with 
the VIRGO simulations having larger simulation boxes and 
lower mass resolutions than the GIF ones. From the VIRGO 
simulations we have analyzed the ACDM model in order to 
test the MJW model in a volume large enough so that the 
effects due to the finite sampling volume may be negligible 



Following the procedure given in section 2.1 we have ob- 
tained the volume- averaged correlation functions up to the 
fourth order from the mass distribution and from the var- 
ious halo samples. Analyses have been performed for two 
different cases. In the first, the higher-order moments are 
calculated at the same time when the dark haloes are iden- 
tified. In the second, haloes are identified at some high red- 
shift while the calculations of the higher-order moments are 
performed for their descendants at a later time. In all cases, 
the redshift at which halo identification is made is denoted 
by zi , while the redshift at which the higher-order moments 
are calculated is denoted by zq. 

In Figures |l|and| we show the third- and fourth- order 
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Model 


n 


n A 


h 




r 


Box Size 


N p 


m p /h~ 1 M 














[Mpc/h] 






GIF-tCDM 


1.0 


0.0 


0.5 


0.60 


0.21 


85 


256 3 


1.0 x 10 10 


GIF-ACDM 


0.3 


0.7 


0.7 


0.90 


0.21 


111 


256 3 


1.4 x 10 10 


VIRGO-ACDM 


0.3 


0.7 


0.7 


0.90 


0.21 


239.5 


256 3 


6.86 x 10 10 



Table 1. Parameters characterizing the simulations used in the paper. Qo and Q A are the density parameters for matter and for 
the cosmological constant, respectively, h is the Hubble parameter, erg is the rms of the density field fluctuations in spheres of radius 
8h~ 1 Mpc, and T is the shape parameter of the power spectrum. Also given are the size of the simulation box, the total number of 
particles and the mass per dark matter particle in a simulation. 



moments from the VIRGO ACDM simulation, together with 
model predictions. For haloes identified and analyzed at a 
given epoch (Figure [j]) both models work remarkably good 
for scales larger than ro [£2 = 1] , which are the validity scales 
of the models, and in particular for haloes with intermediate 
masses. Both models are less accurate for very low mass 
haloes and for very massive haloes. On the other hand, in 
the case of present epoch descendants of haloes identified at 
an earlier epoch, both models work remarkably good in the 
validity scales of the models, as shown in Figure ^. 

In these two figures we also plot the prediction of the 
MJ W model with S3 and S4 given by perturbation theory 
(see Bernardeau 1994). The fact that this prediction also 
matches the simulations results suggests that the moments 
obtained from the VIRGO simulations are not affected sig- 
nificantly by the finite-volume effect and that the MJW 
model is a good approximation to the higher-order moments 
for haloes that are not much smaller than M+ [defined by 
a{K'L) = 1.68]. 

With their high mass resolutions, GIF simulations al- 
low us to test the theoretical models for haloes with mass 
M <C M*. Since the GIF simulations have relatively small 
simulation boxes, the higher-order moments are ex pected to 
be af fected significantly by the finite- volume effect (Colombi 



et al. 1994). However, this effect in each simulation is ex- 
pected to be similar for both the mass distribution and the 
halo distribution. Thus, to test the bias model given in ( |i"o| ) 
and ( pd[ ) by a numerical simulation we should use the value 
of 53 and S4 obtained directly from mass distribution in the 
simulations, because it is the simulated power spectrum (not 
the theoretical one) that is responsible for the clustering in 
the simulations. Figures ^ and ^ show the results obtained 
for the GIF simulations, along with the predictions from the 
models, for present epoch descendants of haloes identified at 
earlier epochs. It can be noticed that the predictions from 
both models, MJW and SMT, are in good agreement with 
the simulations results. 

As we want to test and compare the performance of the 
models for low mass haloes at a given redshift, we have es- 
timated the moments of counts-in-cells, as well as the mo- 
ments predicted by the models, for several populations of 
sub-M* haloes. It is well known that many of the low-mass 
haloes identified using a simple friends-of-friends halo find- 
ing algorithm have positive total energy (i.e., are unbound 
groups) and that these unbound haloes are ex pected to have 



a sig nificant effect on higher-order statistics ( Benson et al 
2001). For this reason we have revised our GIF halo cat- 



least bound particle in the group is removed and the total 
energy is computed again. This process is repeated for each 
halo until the total energy becomes negative or until there 
are remaining less than 10 particles in the group. If the last 
condition is reached, the halo is removed from the catalog. 
We found that ~ 10% of the haloes in the or iginal catalog 
have t o be removed form it, in agreement with 



( [200 1| ) 



Benson et al 



alogues for possible unbound haloes. The total energy for 
each halo in the catalog is computed, if it is positive the 



Figures ^ and ^ show the results for present epoch 
haloes less massive than K'h in the revised GIF halo cat- 
alogues, along with the predictions from the models. For 
haloes with masses much smaller than M*, the SMT ex- 
tension overestimates the skewness and kurtosis, while the 
MJW model gives a better fit to the simulations results. 
Thus, although the SMT extension improves the models for 
the mass function and the second order model of dark haloes, 
it seems to be unable to improve the models for the higher- 
order moments. It is not clear why this happens, it might be 
a fluke or it might also mean that the simple model for the 
higher-order moments, either spherical or ellipsoidal, does 
not work accurately. 

From a comparison between the VIRGO and GIF re- 
sults, it is evident that both the skewness and kurtosis are 
strongly affected by the finite- volume effect. However, if the 
loss of clustering power due to the finite volume is taken 
into account, the model predictions are in good agreement 
with the numerical results, suggesting that the bias rela- 
tions given by ( ji"o| ) and (pd|), with the coefficients given by 
the extended Press-Schechter formalism, are good approxi- 
mations to the skewness and kurtosis of dark haloes in the 
quasi-linear regime. 

To see more clearly the difference between the MJW 
model and the SMT extension, we show in Figure Qthe am- 
plitudes of the halo skewness and kurtosis at a fixed radius 
(R = 10/i _1 Mpc) as a function of the linear bias parameter 
b = bi [see equations @ and @]. The curves correspond 
to the predictions from the models for the present-day de- 
scendants of haloes at three values of Zi (3.0, 1.0 and 0.0). 
From the figure we see that in all cases the values of Sj,h 
are lower than those for the mass unless b is comparable to 
or smaller than 1. This result was obtained in MJW based 
on the spherical model. We see that this is also true even if 
the SMT extension is used, although the amplitudes of Sj,h 
given by the elliptical model are higher than those given by 
the spherical model for a given b. These features in 5*3^ and 
S4 : h have been used in MJW to constrain the bias parameter 
b for galaxies. 
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Figure 1. Skewness S3 and kurtosis S4 of dark haloes with different mass ranges obtained from counts- in-cells analysis (symbols), from 
applying the bias model from MJW (solid line) and its SMT extension (dashed-line). The moments for the mass distribution are shown 
by the dotted line and the moments for the haloes obtained using the moments for the mass from the perturbation theory are shown as 
a dot-long dashed line. The thick ticks on the horizontal axis show the scales where ^2(R) = 1. Results are shown for the VIRGO ACDM 
simulations. The haloes have been identified and analyzed at the times written in the upper-left boxes. The value of M t is also written 
for more information. Each box corresponds to a different range of masses of haloes, as appearing in the labels. The quantities between 
parenthesis correspond to the number of haloes in each sample. 




Figure 2. Skewness S3 and kurtosis S4 for haloes in the VIRGO ACDM simulations for haloes identified at z = 3 and analyzed at 
present time. The lines and symbols correspondence is the same as in Figure hi as well as the notation. 
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Figure 3. Skewness 53 and Kurtosis 54 of dark haloes with different mass ranges obtained from counts-in-cells analysis (symbols) and 
from applying the bias model from MJW (solid line) and its SMT extension (dashed-line). The thick ticks on the horizontal axis show 
the scales where ^2(R) = !• Results are shown for the GIF ACDM model and for haloes identified at z = 1 and analyzed at the present 
time. 




Figure 4. Skewness 53 and Kurtosis 54 of dark haloes with different mass ranges obtained from counts-in-cells analysis (symbols) and 
from applying the bias model from MJW (solid line) and its SMT extension (dashed-line). The thick ticks on the horizontal axis show 
the scales where !;2(R) = 1. Results are shown for the GIF ACDM model and for haloes identified at z = 3 and analyzed at the present 
time. 
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Figure 5. Skewness S3 obtained from counts-in-cells analysis (symbols) and from applying the bias model from MJW (solid line) and 
its SMT extension (dashed-line) of haloes less massive than M*. Each row in the panel corresponds to a different range of halo masses, 
as indicated in the boxes. 



4 DISCUSSION 

From the results shown in Figure ^ we see that for present 
time descendants of haloes already formed at a given redshift 
(z > 0), the values of the skewness and kurtosis depend only 
weakly on the object mass if the bias parameter b is larger 
than or near to one . On the other hand, in the same range 
of b it is clear that the high order moments depend on the 
identification redshift, which is associated with the redshift 
of formation of the objects, with the corresponding values 
increasing as the formation redshift increases. Therefore the 



values of the skewness and kurtosis of old objects, like el- 
liptical galaxies, are expected to be higher than the corre- 
sponding moments of more recently formed objects, such as 
spiral galaxies. This feature can be useful in studying differ- 
ent galaxy populations. 

We have used the models to analyze the predicted val- 
ues of the high order moments for high redshift objects, like 
the Lyman Break Galaxies (LBGs), which are commonly as- 
sumed to for m in the center of the most massive haloes at 
redshift ~ 3 (tvio fc Fukugita |l99d ; |Adelberger et al 



Jing fc Suto 199£ ; Mo et al. 1999 ) . Under this assumption 
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Figure 6. Kurtosis S4 obtained from counts-in-cells analysis (symbols) and from applying the bias model from MJW (solid line) and 
its SMT extension (dashed-line) of haloes less massive than M*. Each row in the panel corresponds to a different range of halo masses, 
as indicated in the boxes. 



and supposing that only a negligible fraction of those haloes 
host a secondary observable galaxy the observed LBGs cor- 
respond to the most massive haloes at z ~ 3. We have es- 
timated the predicted values for the skewness and kurtosis 
at a fixed scale R — 10 ft -1 Mpc of the LBGs (2 = 3) and 
their descendants at a given redshift z. We chose this value 
of R because the mass density in the universe is still in the 
quasi-linear regime and the high order moments of galaxy 
distributions are more difficult to measure on much larger 
scales. 



by equations (fh3"|)-(|l5|), where the S q (q = 3,4) for the mass 
distribution are o btained from linear perturbation theory 
(Bernardeau 1994), and the weighted average needed to get 



Fo r our estimates we have used the coefficients as given 



the effective bk's is derived either by means of the mass func- 
tion from the Press-Schechter formalism, for the predictions 
from the MJW model, or by means of the mass function 
predicted by the SMT model, for the predictions from the 
SMT model . The main parameter for the estimation of the 
bk's for the LBGs corresponds to the observed abundance 
of LBGs, n amely the number density given by (Adelberger 
~ 1998[). This number is Ni bg w 8 x 10" 3 /i J Mpc"' 3 at 
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Figure 7. Predictions from the MJW model (solid lines) and 
its SMT extension (dashed-line) for the skewncss and kurtosis of 
haloes at a radius R = 10 h~ 1 Mpc as a function of the linear 
bias parameter b. each pair of curves shows the results for a given 
Si, where zi = (5i/1.686 — 1) = 0., 1.0, 3.0) from bottom to top. 
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Figure 8. Skcwness and kurtosis at R = 10fo -1 Mpc for the 
LBGs at 2 = 3 and their descendants at later epochs. The curves 
correspond to the ACDM model (dashed lines) and the tCDM 
model (solid lines). The left panel shows the predictions from the 
MJW model and the right panel shows the predictions from the 
SMT extension. 



z ~ 3 for an Einstein-de Sitter universe, and is similar to 
the present abundance of L* galaxies. The corresponding 
number for the ACDM universe is estimated by multiply- 
ing this number by the comoving volume per unit redshift 
at z ~ 3 for an Einstein-de Sitter universe divided by the 
corresponding value for the ACDM universe. 

In Figure |^ we show the values of the skewness, the kur- 
tosis and the linear bias at R = 10 ft -1 Mpc of the LBGs, as 
a function of the redshift, in the ACDM and rCDM mod- 
els. From the curves we see that, although the linear bias 
parameter is quite different in both CDM models, the val- 
ues obtained for the moments are quite similar and so these 
statistics do not provide stringent constraints on cosmolog- 
ical parameters. Note that because LBGs ar e highly biased 
relative to the mass (Adelberger et al. 1998), the skewness 
and kurtosis parameters of their haloes and descendants are 
significantly lower than those for the mass. Thus, observa- 
tions on these quantities may give additional evidence that 
these objects are highly biased. 

Finally in Figure (9) we show the higher-order moments 
from the spatial distribution of model galaxies in the GIF 
ACDM simulations at redshifts {z = 0, 1 and 3), along with 
the corresponding quantities of the distribution of haloes 
hosting the galaxies and with the same quantities of the 
distribution of haloes weighted by the number of galaxies 
hosted by each halo. The catalogues are limited to model 
galaxies with stellar masses greater than ~ 2x 10 10 h~ 1 Mq 
(i.e., haloes more massive than 10 12 /i -1 M @ ). For fur- 
ther details about these catalogues and t he galaxy forma- 



t ion m odels used in their construction see [Kauffmann et al 
( [1999D 



From the plots it can be seen that the skewness and 
kurtosis of model galaxies at redshift 1 and 3 is similar to 
the corresponding quantities of the dark matter halo popu- 
lation hosting the model galaxies (Mh > 10 12 h^ 1 Mq). The 
difference between the moments of the model galaxies and 
dark matter haloes are still appreciable since there are haloes 
hosting more than one galaxy and this have strong effects 
on clustering quantities like the higher-order moments. Nev- 
ertheless, the similarities shown in the plots suggest that, 
indeed, one may assume that the distribution of galaxies 
may be approximated by the distribution of dark matter 
haloes. Indeed from the plot it can be seen that the mo- 
ments of the distribution of model galaxies are very similar 
to the moments of the distribution of dark matter haloes 
weighted by the number of galaxies hosted by each halo. At 
small scales the moments of the weighted halo distribution 
are larger than the moments of the model galaxies distribu- 
tion because we neglect any structure information in a halo 
containing several model galaxies. 

Thus, one might be able to describe the higher-order mo- 
ments of galaxies by combining the theoretical models for 
the higher-order moments of dark matter haloes with mod- 
els for the number, stellar mass and position of individual 
galaxies that can be hosted by a single dark matter halo 
with mass M at a given epoch. Furthermore, many of the 
semi-analytic model galaxies at redshift 3 would be included 
in current LBGs samples and thus are close to a subset of 
the LBGs population, which is highly biased respect to the 
mass. From the figure we see that S3 ~ 2 and 54 ~ 5 at 
R = 10 Mpc/h, which is nearly consistent with the model 
predictions shown in Figure H 
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Figure 9. Skewness S3 and kurtosis S4 of the spatial distribution 
of galaxies obtained from the GIF ACDM simulations using semi- 
analytical models of galaxy formation (filled circles). S3 and S4 
of the distribution of haloes hosting the galaxies is also plotted 
(open circles). For comparison, the respective quantities of the 
distribution of haloes weighted by the number of galaxies they 
host is included (crosses connected through a line). 



5 SUMMARY 

We have tested the MJW model and its SMT extension using 
two sets of high-resolution N-body simulations with differ- 
ent simulation boxes and mass resolution. From the set with 
the very large simulation box, which allows us to control the 
finite volume effect, we found that the models work remark- 
ably good for CDM universes. The good performance of the 
models when the moments from the mass distribution are 
estimated using the linear perturbation theory, shows that 
the moments from this set (VIRGO) are practically unaf- 
fected by the finite volume effect. The other set of simula- 
tions, having much higher mass resolution, has been used to 
test the model for low-mass haloes, showing that the model 
based on spherical collapse (MJW) works better than the 
model based on ellipsoidal collapse (SMT) in describing the 
higher-order moments of haloes less massive than M* . For 
massive haloes both models work remarkably good. 

We use the theoretical models to predict the higher- 
order moments at a fixed scale of the Lyman break galax- 
ies observed at z — 3 and their descendants at lower red- 
shifts. We found that, although the linear bias parameter b 
depends strongly on the cosmology adopted, the values of 
the higher-order moments are practically the same in both 
CDM models and therefore the higher-order moments from 
the spatial distribution of these objects cannot be used to 
constrain cosmological parameters. 
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